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Abstract 



We study the use of polarization asymmetries at a muon colhder to determine the 
CP-even and CP-odd couphngs of a Higgs boson to fi~^fi~ . We determine achiev- 
able accuracy as a function of beam polarization and luminosity. The appropriate 
techniques for dealing with the polarization precession are outlined. Strategies 
especially appropriate for a two-Higgs-doublet model (including the MSSM) are 
given. Our general conclusion is that polarization will be very useful, especially if 
the proton source is such that full luminosity in the storage ring can be retained 
even after imposing cuts on the originally accepted muons necessary for P > 0.4 
for each beam. 
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1 Introduction 



Although the origin of mass has still not been estabhshed, it is widely expected that 
electroweak symmetry breaking is driven by elementary scalar dynamics, leading 
to one or more physical Higgs bosons [jl| . It is also very possible that CP violation 
arises either partially or entirely as result of CP violation in the Higgs sector [Q]. 
Even if not, CP violation in other sectors of the theory can induce CP violation 
in the Higgs sector at the loop level 0. Thus, we anticipate that the direct 
determination of the CP nature of each observed Higgs boson could be crucial to 
unraveling the nature of the full theory. 

Even if the minimal one-doublet Standard Model turns out to be nature's 
choice, we will certainly want to know that the single observed Higgs boson is 
entirely CP-even in nature. Alternatively, if electroweak symmetry breaking turns 
out to be driven by technicolor-like dynamics, it would be highly desirable to be 
able to check that the narrow, light pseudo-Nambu-Goldstone bosons (PGB's), 
that often arise in such a theory and are rather Higgs-like in many respects (see, 
for example, 0), do indeed have CP-odd coupling to fermions. 

The means for determining the CP nature of an observed narrow resonance are 
limited. If it is determined that its WW or ZZ coupling is substantial (small), 
then we will know that it has a substantial (small) CP-even component. But, even 
in the general two-Higgs-doublet model the relation between this coupling and 
the CP-even and CP-odd fermionic couplings (denoted a and ^675, respectively, in 
/(a + i&75)/Q) is model-dependent. However, polarization correlations can allow 
one to extract the ratio 6/a in a model independent manner. For a light resonance, 
one possibility is to employ the r+r~ decay of the resonance. If the tt mode is 
open, then it too can be employed. These final-state possibilities were examined 
for a muon collider in p, p. One can also look for certain characteristic angular 
distributions in associated production, 66-|-resonance or resonance, that are 
sensitive to b/a [0. However, the most elegant approach to determining the CP 
nature of a neutral resonance is to employ initial state polarization asymmetries. 
This is possible only for 77 or fi'^fi~ [§], ^ production of the resonance. But, the 
77 coupling of a neutral resonance is the result of either loop(s) (in the Higgs case) 
or an anomaly (in the PGB case), and is not a direct measure of the elementary 
fermionic couplings. This leaves /i^yU^ collisions, which are also the only way 
in which we will probe 2nd generation fermionic couplings in models where the 
strength of the fermionic coupling is proportional to the fermion mass (implying 
that a resonance with mass > 2mr will always decay to 3rd generation fermions 
and any CP-odd coupling to two photons will be dominated by the top quark loop). 

In the limit of = ^1 — 4m^/m|. 1, the cross section for production of a 

"^Phases can always be chosen so that, for any given /, a and b are reaL 
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resonance, R, with fi{a + ibj5)fi coupling to the muon takes the form 



wsK) =i^s[i + PIPE + PfPE - ^ 

= ifj [l + P+PE + P+Pf cos(2i + C)] , (1) 

where 6 = tan~^ - and Pt (Pl) is the degree of transverse (longitudinal) polariza- 
tionQ of each of the beams defined as P = f^ — f~, with being the fraction 
of muons with spin in the dominant direction and /~ the fraction of muons with 
opposite spin. Here, ( is the angle of the /i+ transverse polarization relative to 
that of the fi~ as measured using the the direction of the /i~'s momentum as the 
z axis. The S subscript denotes 'signal'. Both the cosC and sin^ dependences 
allow significant sensitivity to the ratio of interest, b/a, even though only the sm( 
term is truly CP-violating. The value of a^, which results from convoluting the 
resonance shape with a Gaussian distribution in ^/s, depends upon the model, the 
detector, the final state mode and the machine parameters. Ignoring final state 
efficiencies and acceptance cuts, and integrating over final state phase space, one 
finds a result that can be approximated byQ 

-0 47rr(P -> fi+fi-)BR{R -> F) 

'^s = f ; ^ 2^ 1/2 ' (2) 



where is the Gaussian resolution in y/s. 

It is expected that the muon collider will first be operated with the relatively 
small natural polarization of P ~ 0.2 for the /i"*" and fi~ bunches, since this will 
allow maximal machine luminosity and lead to the largest number of Higgs events. 
The bunch polarizations will be oriented in the plane of the final storage ring so 
that their precession as the muon bunches circulate (for roughly 1000 turns) can 
be used to precisely determine the central energies of the muon bunches and their 
Gaussian energy spread JlT], |12|. (One measures the oscillations of the energies of 
the electrons from the decaying muons, which oscillations depend very sensitively 
upon both the muon energy and the Gaussian energy spread.) From Eq. ([^) and 
the fact that |P^P^| and \P^P^\ are both < 0.04, it is clear that these obser- 
vations will, to an excellent approximation, provide a measurement of a^. (By 
choosing appropriate relative phases between the /i"*" and /x~ polarization angles, 
the polarization dependent terms in Eq. (||) will average to zero — see the next 
section.) Thus, a very accurate measurement of in several channels F can be 
performed and a reasonably accurate measurement of F^* by a three-point scan 

"^Note that the cross section for longitudinally polarized muons depends only on + and 
cannot be used to extract information about the CP nature of the resonance's couplings to muons. 

"^This form, fr om Ref. has the correct asymptotic limits for Tj^/a^ and Tji/a ^ 
oo and is always within 18% of the exact result. 



2 



will be possible. A model-independent determination of + }? then becomes 
possible from Eq. (H) by computing r(i? — >• using the measured values of 

0^, r*^*, and cr^ and a determination of BR{R F). To obtain BR{R F) 
requires using the missing mass technique in the ZR final state (in either e~^e~ 
or fi~^fJ^~ collisions) to measure a{ZR ZX) {X = anything) and the measure- 
ment of a{ZR ZF) to compute BR{R ^ F) = a{ZR ZF)/a{ZR ZX). 
Once the basic Higgs observations have been performed so that the Higgs width, 
branching ratios and -|- 6^ are well-determined, the next task will be to perform a 
model-independent measurement of b/a. This require maximizing the influence of 
the polarization-dependent terms in Eq. (|l|). Determining the best procedure for 
doing so and estimating the accuracy with which this measurement can be carried 
out is the main goal of this paper. 

Any given final state F will have a significant background. In general, the 
cross section for the background is nearly independent of a^. If we integrate 
over final state configurations then the background is independent of P^t and 
and (. In order to roughly understand the level of sensitivity to b/a that can be 
achieved, let us for the moment imagine that we can set = Pjf = and choose 
Pjt = P^ = Pj.. We then denote the integrated background cross section by a%. 
We imagine isolating and -^f^, respectively, via the asymmetries 

-4,. !^ -'i =Pj^ = P|cos2.. (3) 

(TsiC = 0) + asiC = vr) + ¥ 

_ as{C = rr/2)-as{C = -n/2) _ 2ab _ , 

as{C, = 11/2) + as{C, = -11/2) + ¥ 

Assuming that (a^ + 6^) must be measured at the same time as the asymmetries, 
the error for the measurement of either of the asymmetries Ai or An is given by 

^^^^ " lMp 

where L is the integrated collision luminosity, assumed distributed equally between 
the two ( measurements and we have temporarily taken Pp to be the same for all 
collisions (which is not actually the will be discussed shortly). If (a^ + 6^) 

has already been very accurately determined using initial resonance production 
measurements and the technique outlined below Eq. (Q), then the expression for 
[SA]"^ takes the form 



i^-^l = n=oi^ (6) 



In either case, if Py < 0.5, and since A"^ < P^, the A^ in the numerator of [SA]"^ 
can be neglected, implying that is proportional to P^L. 

From this discussion, it is apparent that the ideal situation would be to arrange 
for the muons to be entirely transversely polarized at the interaction point (IP) 
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and to be able to adjust the angle between the /i~ and polarizations to be 
fixed at one of the four values ( = 0, 7r/2, tt, 37r/2 for each collision. However, 
this is not possible if we are interested in a very narrow resonance, such as a 
light SM-hke Higgs boson which will have a width of just a few MeV. The reasons 
follow. First, it is important to note that without intervention any horizontal (i.e 
perpendicular to the magnetic field of the storage ring) polarization will precess 
about the magnetic field and therefore rotate relative to the momentum direction. 
Consequently, the amounts of transverse and longitudinal polarization at the time 
the bunch passes the IP will oscillate. This will be described in more detail in the 
following section. The only way to overcome this oscillation would be to have a 
section of the storage ring devoted to compensating for this spin precession on a 
turn-by-turn basis. However, this conflicts with the requirement that we be able 
to measure the central energy of each muon bunch circulating in the storage ring 
to 1 part in 10^ (and the beam energy spread of the bunch to better than 1 part in 
10^) using measurements of the oscillations of the energy of the secondary electrons 
from the muon decays as the bunches circulate and their spins precess. A large 
number of turns in which the precession is allowed to occur without compensation 
is required. The above precisions are those needed for a light Higgs resonance 
with a SM-hke width (of a few MeV) in order to be certain of remaining rather 
precisely centered on the resonance peak and knowing exactly how that resonance 
peak is being sampled. For a resonance with width of a few hundred MeV or 
larger, knowledge of the beam energies to 1 part in 10"^ would be adequate. Since 
magnetic fields would be known and be stable at this level, direct measurement 
of the bunch energies would not be required and one could consider turn-by-turn 
spin compensation so as to achieve the ideal transverse configurations at the IP for 
each collision. Here, we will assume that the storage ring will not initially be built 
with the extra magnetic components required for such compensation. 

Thus, regardless of whether the Higgs or other resonance is broad or narrow, it 
is necessary to determine the effects of spin precession on the accuracy with which 
cos 25 and sin 25 can be determined. In what follows, we develop a procedure 
whereby the luminosity needed to achieve a given accuracy is only about 50% 
larger when the spins are allowed to precess than if the spins could be taken to be 
purely transverse for each collision. 

2 Polarization Precession at a Muon Collider 

Most muon collider designs are such that there are two bunches, each, of //+'s and 
//~'s circulating in opposite directions in the storage ring. Typically, these bunches 
will be stored for about 1000 turns. Once the bunches enter the storage ring, any 
polarization in the horizontal plane will precess in the vertical magnetic field of the 
storage ring. Further, as noted in the introduction, for the case of a very narrow 
resonance this precession is needed to measure the bunch energies and Gaussian 
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Figure 1: We plot P^{Nt) I\Ph\ at the interaction point as a function of the 
number of times, Nt-, that the /i~ beam passes the IP, assuming the muon coUider 
is operating at a total center of mass energy of a/s = 110 GeV and that the /i~ 
enters the storage ring with longitudinal polarization p}j/\Ph\ = 1- This plot is 
from M . 



spreads to high accuracy, implying that the polarizations of the bunches should not 
be manipulated once they are stored in the ring. The impact of such precession 
on extracting physics has not been carefully examined to date. Thus, we provide 
in this section a fairly detailed explanation of the considerations and procedures 
that must be employed. Even though we shall focus on the case of a narrow spin-0 
resonance, the general features of our discussion will have wider applicability. 

The typical configuration will be such that a given bunch enters the ring with a 
component Py {Pr) of polarization vertical (horizontal) with respect to the plane 
of the storage ring. [Py and Pr are defined in the muon rest frame; see Eqs. (0) 
and (H) below.] P^ and will rotate in the same directions as the /i~ and /x"*" 
themselves (that is in opposite directions). The rate of rotation of the Prs as 
viewed from the laboratory frame is somewhat different than the rate at which 
the bunches themselves rotate. The mismatch means that if, for instance, P^ 
were longitudinal at the time the /i~ bunch first enters the storage ring, it will not 
remain so but rather it will precess into the transverse direction and then back to 
the opposite longitudinal direction, and so forth. 
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To be precise, we assume that the storage ring's magnetic field points in the 
—y direction: B = —By. The n~ is then rotating clockwise about the y axis. At 
any given moment, we define z = p^-. Then, x points radially outward. The angle 
9~ is defined as the angle by which one must rotate (in the /i~ rest frame) about 
the y axis in order to get from z to ; 9~ is an 'over-rotation' angle in that it is 
the additional angle of rotation of the spin as compared to the angle of rotation 
of p^- . After boosting to the laboratory frame, the complete four-component spin 
vector for the /i~ is then: 



Spi Ph 



7(/5,£)cosr - (0,J)sinr +P^(0,y). (7) 



The standard result of Ref. |T3[, assuming that the /i~ beam enters the storage 
ring with P^ = p^-, is 6~{Nt) = uj{Nt — 1/2), where Nt is the number of turns 
during storage, counted starting with A^^ = 1 the first time the bunch passes the 
IP {Nt = 1/2 at bunch insertion), and u = 271'^^^^^, with 7 = E/m^ and ^^^^ = 
1.165924 X 10"^ One finds that = vr for a beam energy of 90.6223 GeV/2 (i.e. 
close to m^/2), implying that an initially longitudinal /i~ horizontal polarization 
becomes transverse after the /i~ travels half way around the ring to the interaction 
point. However, for the somewhat higher energies of interest for a Higgs factory, 
the degree of horizontally transverse polarization (denoted by IS) oscillates with 
Nt according to Ph{Nt) /\Ph\ = sin 9~{Nt) as illustrated in Fig. where we have 
chosen a convention in which P^ is positive when it points towards the center of 
the storage ring. We now discuss how to take this oscillation into account. 

First, we note that similar results apply for the /x"'". The will be traveling 
in a counter-clockwise direction about the y axis; at the interaction point 1/2 way 
around the ring, p^+ = —z. Once again, we can define an over-rotation angle 9^ 
(in the /i+ rest frame), in terms of which the /i^'s spin vector (in the laboratory 
frame) at the interaction point is written as 



7(/5,-2)cose+ - (0,£)sin^+ +P^{0,y). (8) 



If we start with P^ = p^+ at the time of insertion, then 9~^{Nt) = uj{Nt — 1/2), 
just as for 9~{Nt)- More generally, we can insert the /i+ beam with any initial 
angle for P^ that we desire. 

We assume that each time a Higgs event is observed we can compute, or will 
have measured (prior to the IP, and then extrapolated to the IP), the transverse 
polarizations of the bunches. That is, we will know 9^ and 9~^ for each interaction. 
In fact, 9^ will be completely correlated with 9~ , the correlation being determined 
by the initial spin configuration with which the bunches are injected into the storage 
ring. We now give the expression for the cross section as a function of 9~ and 9~^, 
defining c_ = cos6'~ etc., 

""'^^^^ ^ = (1 + P^PhC+c^) + cos25{P^Py + P^jPhS+s.) 

+sin25(P^P+s_ - P^P^s+) . (9) 
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Note that we obtain Eq. (|l|) using the obvious replacements: PhC+ = Pl, Pr^- = 
PE, P^Py + P^PhS+s_ = P^Pt cosC, and PffS-P^ - PrS+Pv = -PtPt sinC- 
We also note that if we choose Py = Py = and insert the bunches so that 
6^ = 9^ + 7i/2, then all polarization dependent terms in Eq. (P) average to zero 
after many turns. This is the configuration that would normally be employed for 
the initial Higgs resonance scans so that knowledge of b/a would not be needed in 
order to properly interpret these measurements. 

However, to determine b/a we will wish to employ one or more different polar- 
ization configurations and retain maximum information by binning Higgs events 
according to {6^ ,6^{6^)). After integrating over {6,(f)) configurations in an // 
final state, one finds (in the limit of rrif/^/s 0) 

anie-, r ) (X ejej(l - P+P^c+c_)H; 

+ 2e/e^c/ [c^((l - P+PhC+c^) + d^^PnC- - P^c+)\ H^3ft(Hz) 
+ (4 + 4) [(cj + dl){l - P^PhC+c_) + 2cMPhC^ - P+c+)l |H 



|2 
Z\ ; 

(10) 



where the 7, Z propagators and couplings are H^ and Hz and icf and i'^f^^Cf+dfy^), 
respectively; see the Appendix for more details. 

In order to give a simplified discussion, let us assume that a| (and, thence, 
a^ + b'^) has been precisely determined by first resonance production measurements 
that do not focus on determining the Higgs CP properties. As we have described 
in the introduction, this is very likely to be the case. Then, remembering that 

can be considered to be a function of 6~ for any given choice of storage ring 
insertion configuration, we can write the signal cross section as 



^sin = /o''(r) +cos25/f (r) +sin25/f (r) , (11) 



where 



f^{e-)=a'si^ + P^PHC+c.), 
fciO-) = ^s^P^Pv + P^PhS+s_) , 

/f r ) = ^s{PhPvS- - P^PyS-,) (12) 

depend upon the configuration, C, chosen for P^, P^, Py , Py ,9^(9^). (Recall 
that 9~^{9~) is determined by the choice made for the relative angle between the 
polarization of the /i^ as compared to that of the fi~ at the time of bunch insertion.) 
The assumption that + 6^ has been precisely determined already corresponds 
to assuming that /o"(^~) is completely known. Then, writing Ti'-^{9^) = o'g{9^) + 
a§(^~), the Ax^ difference between two different Higgs models with the same 
+ 6^ but different values of S will be given by 

Ax' = Y.Lc [^iMg + 2A,AsMg + A^m£] , (13) 
c 
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where Lc is the luminosity devoted to configuration C, 



m: 



c_ fd9- if^je-)] 

J 2n 



m: 



c 
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M. 



c 



de- [/f (r)] 

27r T.c{e- 



(14) 

and we have defined Ac = A cos 25 and A^ = A sin 25 (the differences in these two 
quantities between the two models) .0 The 6~ integral is over the value of 6~ (and 
the correlated 6^ (9^) value) at the interaction point. Here, we are approximating 
the sum over the discrete 9^ values that arise during the course of 1000 turns 
of the bunches by an integral over 6^ . This is an excellent approximation unless 
-v/i/90.62 GeV (where 90.62 GeV is the i/i value such that 9~ is constant) is a 
ratio of small (compared to 1000) integers. 

We will demonstrate below that the configurations can be chosen sufficiently 
cleverly that one can neglect the 9~ , and, indeed, the entire configuration depen- 
dence of T,'-^{9~) and approximate S'"(6'~) ~ + 0%, the sum of the unpolarized 
cross sections. Under these circumstances, the Ax^ for discriminating between two 
different models with the same + 6^ value (i.e. same total unpolarized rate) will 
then be given by 



2 

C ' 



B C 



where 
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A,{P+Py + P+PhS+S_) + A.,{P^P+S_ - P^PyS_ 



(15) 



(16) 



c 



is a measure of our 'sensitivity' to cos 25 and sin 25. The averaging is over roughly 
1000 turns of the bunches in the storage ring. To an excellent approximation, this 
average depends only on the relative values of 9^^ and 9~ at the time of bunch 
insertion. We will now give results for for (a) 9^ = 9^ (b) 9^ = 9^ + ti (c) 
9- = 9+ + Ti/2. The resuh for (d) 9- = 9+ + ?>ti/2 is the same as for (c). We find 



5! = Al 



fi2 



S": = Al 



-{P^Pnf + P^PhP^Pv + {P^P^-'' 
3 



liPiiPnf 
1 
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{PiIPH? + {P^P^ 



Pr Pr Py Py + (-^y Py 
1 



y ) 

-^2 



yf 



2^ 



P^Py 



PrP^ 



1 

2' 



+ ^ A^ \P^Py + P^P + 



+ 2^^ 



{P^Pyf + {PhP^Y 



(17) 



We note that M^s = when the 9~ dependence of Tp is neglected. 

As noted earlier in the previous section, the ideal situation would be to be able 
to choose four basic relative orientations of the /i~ and /x"*" transverse polarizations 

"^We note that in the hmit of just two bins corresponding to sin^ = ±1 or cos^ — ±1, A^^ = 1 
corresponds to errors for sin 25 or cos 25, respectively, as given in Eq. (0). 
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with respect to each other: C = 0; ^/2, tt, 37r/2. Constant C = or ^ = vr could be 
accomphshed by setting = = and Py = Py or Py = —Py, respectively. 
However, as noted in the previous section, some degree of horizontal precessing 
polarization is necessary if we are to be able to measure the energies of the muon 
bunches with the accuracy of 1 part in 10® needed for a very narrow Higgs reso- 



nance. It is estimated that Pj^ = Pjj ~ 0.05 — 0.1 IjT^ is required for the energy 
measurement. (We will adopt the optimistic choice of 0.05 in the remainder of this 
paper.) The rest can be placed in the vertical directions. 

The spin precessions make it impossible to maintain C, = vr/2 or ( = 37r/2 
as the bunches circulate. The simplest thing that one can do is to inject, say, 
the /i~ bunches with purely horizontal polarization, and maximize the vertical 
polarization for the /i"*" bunches subject to the requirement that P^ = 0.05. Sen- 
sitivity to both the magnitude and sign of sin 25 arises from the variation of the 
sm26P^S-Py term in Eq. (|16[), which at various extremes samples ( = 7t/2 and 
( = 37r/2. (We again emphasize the importance of either calculating, starting from 
the insertion configuration, or measuring, from the decay spectrum, the precession 
angles associated with each observed Higgs event.) 

Very specifically, we thus consider the following configurations, keeping in mind 
that the net polarization P for the fi~^ and fi~ bunches will be essentially the same. 

I: To approximate the C = configuration, we choose P^ = P^j = Ph = 0.05, 
9~ = 9~^, Py = Py = \J P"^ — Pjj . The sensitivity is then given by 

§} =Sl = Al [3P^/8 + PUP' - Pi) + {P' - Phf] . (18) 

II: To approximate the = vr configuration, we choose P^ = Pjj = Ph = 0.05, 



9 = 9^ + n, Py = —Py = —^P"^ — Pff. The sensitivity is then given again 
byEq. (|8D: = 5|. 

We now justify, for the case of C = J, II, the approximation of taking S*^ ~ 
a^-|-a^, employed in obtaining Eq. ( |T5|) from Eqs. ([13|) and (0). We note that for 



both configuration I and configuration II the background, see Eq. (|TOD, will depend 
only very weakly on {9~ , 9^) simply because P^ and P^ are both small. Similarly, 
the P^PJjCj^C- term in f^ij)') and the P^P^s+s^ term in will be very 

small. Further, f^{9~) can be approximately neglected because both P^ and P^ 
are small. Finally, we will sum over configurations I and II (with equal luminosity 
weighting) in order to determine cos 25. This has the important consequence that 
(using \Py\ — P, see above) the leading term in Eq. (|TB|) is obtained from 
Eqs. (|T3|) and ( [T^ ) via the structure 

EnrC p4r— 0l2 [ ]_ I ^ ] 
c^jjj ^ \a% + a%{l + P^cos25) a% + a%{l-P^cos25)) 



2(^B+^' 
(a0j + a0)2-P4cos2 25(^0^2 
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so that even for P as large as 0.5 or so (as we shall later consider) the correction 
term in the denominator can be neglected.0 Thus, the approximation of taking 
E*^ ~ 0^ + is quite accurate in this case. 

Ill: To emphasize the ( = 11/2 and ( = 37r/2 configurations over many turns of 
the bunches, we choose = P {Py = 0), P^ = Ph = 0.05 and Py = 
^JP^ - P2 . In addition, we choose e~ = 6+ + 7r/2 (or 6' = 9+ + 3n/2) so as 
to minimize dependence on Ac. The sensitivity is then given by 

S]n = Sl = l^lP'Pl + l^lP'iP' - Pi) ■ (20) 

Note that if we choose either 9~ = 6^ or 6~ = 9^ + vr, then the (undesired) 
sensitivity to Ac would be increased. For Ph < 0.05 and P > 0.2 (as will be 
the case), the A^ term can be dropped in Eq. (^Oj). 

In order to justify neglecting the dependence of S*^ on 6~ in case III, we note that 
the terms of concern in aB{0~ ,9~^), Eq. (p!0D , are those proportional to {PjjC- — 
PhC+) ~ Pc-. The term of concern in as{0~ , 9^) is that proportional to fs{d~) ~ 
a^P^s^. We will return to this latter term in a moment. First, we note that 
the Pc_ term in aB{O~,0~^) can be approximately eliminated without affecting 
as{6~ , 0^) by binning together into a single c° bin events with the same sign of s_ 
but with c_ = c° and c_ = — c° . (Even after 1000 turns, this is only an approxima- 
tion since the number of turns for which c_ ~ c° can differ significantly from the 
number of turns for which c_ ~ — c° , depending upon the Higgs boson mass and 
the number of c_ bins employed.) Concerning the as{0~ , 0~^) term proportional to 
P^S-, we note that in computing M^g, Eq. ([I^) , the leading term in the numerator 
is even in s_ while the P^s_ correction term in the as{0~ , 9^) contribution to TP 
is odd. This means that we will get a cancellation analogous to that in Eq. (p!9|) 
so that these corrections to our approximation will be of relative order P^ and can 
be neglected. 

It is useful to note how these results compare to the ideal where we imagine 
that C can be held fixed. For ^ = 0, vr, our sensitivity S"^ would be A^P"^ while for 
C = 7r/2,37r/2 it would be A^P^. For Ph = 0.05 we suffer a loss of about 0.469 
(0.939) for A, (Ac) for P = 0.2 and 0.495 (0.990) for P = 0.5. In the remainder of 
this paper, we approximate these Ph = 0.05 results by assuming no loss in for 
Ac and a factor of 1/2 loss in S"^ for A^. This latter factor must be overcome by 
increased luminosity to obtain the same statistical accuracy as in the ideal case. 
To equalize sensitivity to A<, and Ac, we will assume that we accumulate twice 
as much luminosity for the A^ configurations as for the Ac configurations. Thus, 
for total integrated luminosity L we accumulate L/6 in configuration I, L/6 in 
configuration II, and 2L/3 in configuration III. Of the latter 2L/3, as the spins 

"^Note the analogy to Eq. (||) that is apparent when we recaU that \SA\'^ corresponds to Ax^ — 1 
and that in the present case is P"^ cos^ 26. 



10 



precess L/3 will be accumulated in configurations with sin^ < and L/3 with 
sm( > 0. The sensitivity achieved is then equivalent to probing the four fixed ( 
configurations {( = 0, 7r/2, vr, 37r/2) with L/6 each. The results of the following 
sections will sometimes be phrased in this latter language. 

3 Maximizing Sensitivity to CP Violation 

The number and density of muons in each bunch are limited by space-charge and 
muon beam power considerations. Existing designs have proton source intensity 
such that these limits are saturated when momentum cuts on the initially accepted 
muons coming from the target are such that the muon beams have polarization of 
P ~ 0.2 (or less). To date, there has been no strong reason for designing the proton 
source so as to have full luminosity (fully saturated bunches) for larger P values. 
This is because the typical design for the storage ring is such that the polarization 
is longitudinal at the interaction point half way around the ring, motivated by the 
fact that (for longitudinal P) the signal rate is proportional to 1 + P^, while 7*, Z* 
induced (e.g. bb) backgrounds are proportional to 1 — P^. In this configuration, the 
improvement of S/ y/S + B with increasing P (assuming the integrated luminosity 
can be kept fixed) is rather slow, rising by < 20% from P = 0.2 to P = 0.4. Thus, 
although higher (longitudinal) P would have some advantage if the proton intensity 
were such that the two muon bunches were both fully saturated at P > 0.4 rather 
than at P ~ 0.2, it is usually accepted that this advantage is not sufficient to 
justify the expense associated with building a more intense proton source. 

However, if determining the CP nature of the resonance's couplings to 
is the goal, the (approximate) proportionality of the S"^ sensitivity measures (see 
previous section) to P^ implies that retaining maximal luminosity at P > 0.4 has 
much greater advantages. The extra proton source intensity required to achieve 
this is determined by the fraction, fsuwiP), of muons that survive after imposing 
the momentum cuts required to achieve a certain P (prior to inserting the muon 
bunches into the storage ring). The variation of /surv as a function of P has 
been estimated by many different groups (see, e.g., [|11|, |12[)- Here, we adopt the 
post-cooling result given in [jl^] (Fig. 21) which is approximately described by 
/surv(-P) = 1-45 — 2.51P for 0.2 < P < 0.6. The number of muons in each bunch, 
Nb, is proportional to the proton source intensity, £ps, times /surv Using the 
specified form, one finds, for example, that /surv(0.45) : /surv(0.39) : /surv(0.2) ~ 
1/3:1/2:1, implying that £ps at P = 0.45 (0.39) must be 3 (2) times that for 
P = 0.2 to maintain full bunches. 

If the machine is constructed with Cps only large enough to saturate the bunches 
at P = 0.2, then one finds that /surv is such that S/ ^/ S + B declines if one increases 
P above 0.2 (and arranges the polarization to be longitudinal at the interaction 
point). However, maximum S is rather achieved by increasing P to the point 
(P ~ 0.39) where each muon bunch has one-half of the P = 0.2 number, N^, and 
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then merging the two bunches into one bunch (prior to insertion into the storage 
ring). Bunch merging effectively doubles the storage ring luminosity at this P 
compared to what would result without bunch merging. (Without bunch merging, 
one would have luminosity proportional to 2 x [Nb/2y whereas with bunch merging 
the luminosity would be proportional to 1 x N^.) As a result, the collider luminosity 
at P = 0.39 is only a factor of 2 lower than at P = 0.2 and P^^/L is roughly a 
factor of 2.7 larger. (At the same time, S/\/S + B has declined by about a factor 
of 1.4.) If CP studies are the goal, one will choose the P = 0.39 option with bunch 
merging. If Cps larger than that required to saturate the muon bunches at P = 0.2 
is available, bunch merging may or may not be desirable for CP studies. For Cps 
such that bunch saturation is achieved for P < 0.42 (> 0.42), the sensitivity S for 
CP studies will (will not) benefit by increasing P still further to the point where the 
bunches can be merged. If the bunch-merging option is appropriate, the integrated 
collider luminosity L will be one-half the (no merging) value at P = 0.2. If bunch 
saturation is the better option, L will be exactly the same as the P = 0.2 value. 
In either case, increasing P beyond the bunch saturation or bunch merging point 
(whichever is more optimal) causes S to decline. 

To be more explicit, we define the relative sensitivity Rs = S/Sq, where the de- 
nominator is that achieved using P = 0.2 and £pg such that the two bunches are sat- 
urated for P = 0.2. Let us also define the relative source intensity / = Cps/ Cp^. In 
Fig. 0, we plot Rs obtained either by bunch merging or bunch saturation (whichever 
gives the larger Rg) as a function of /. Also shown is the corresponding polariza- 
tion being employed as a function of /. The corresponding functional forms are: 
for bunch merging, P™(/) = cl — ^ with R^{I) = ("^75^) 2~^/^; for simply filling 

two bunches, P^(/) = a - | with P;^(/) = {^^f- Here, a ~ 0.577 and b ^ 0.377 
for the form of fsuw given earlier. Below (above) / = 2.38 we employ P"^ and Rg 
(pf and P5). The results for Rg show clearly that if determining the CP nature 
of a resonance is the goal, then one gains by optimizing P and having as large a 
proton source intensity as is feasible. 

Of course, the absolute accuracy with which a measurement of b/a can be made 
is dependent upon the resonance model. In the following section, we will consider 
several Higgs boson examples and demonstrate that by using optimal polarization, 
as described above, a very meaningful measurement of b/a can typically be per- 
formed. Overall, we will conclude that the increase in our ability to determine 
the CP nature of the muonic couplings of a resonance would provide substantial 
motivation for spare proton source intensity. 

Before ending this section, we note that there is a very interesting possibility 
for increasing the number of muons with high polarization retained after making 
the necessary selection cuts. In particular 0, |l^, [1^, if the accelerating gradient 
of the phase-rotation device that immediately follows the pion capture solenoid is 
sufficiently high (~ 4 — 5 MV/m), the correlation between the muon arrival time 
and average muon polarization might be significantly enhanced, resulting in a more 
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Figure 2: We plot the maximum achievable relative sensitivity of the CP deter- 
mination as a function of the relative proton source intensity /, both being relative 
to values corresponding to two full //"^ and iJ,~ bunches (each) at P = 0.2. At 
/ = 2.38, one switches from bunch merging to bunch filling. Also shown is the 
polarization being employed for a given I. 
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effective selection of polarized muons, and a higher final polarization for the same 
luminosity. Such high-gradient phase-rotation designs are being actively considered 
for the neutrino beam facility version of the muon accelerator and storage ring [ITB . 



4 Two test cases: a light Higgs with SM-hke bb 
rate and a degenerate H^-A^ MSSM pair 

Consider first the case of a Higgs boson with SM-like bb final state rate and total 
width. For example, it might be that a Higgs boson is detected and appears to 
have SM-like WW/ZZ couplings and branching ratios to accessible final states, in 
which case we will wish to determine if its fermionic couplings are indeed CP-even. 
Here, we assess the accuracy with which this verification can be performed at a 
muon collider. As detailed in earlier studies P, 0], the result of Eq. (H) and the 
very small width of a light SM-like Higgs boson means that the muon collider must 
be operated with the smallest possible beam energy resolution, even if this results 
in substantial luminosity sacrifice. For currently understood designs, the best that 
can be achieved is Ai^bcamZ-^beam = R with R of order R = 3x 1O~^.0 For such an 
R value, the yearly integrated luminosity is anticipated to be of order L = 0.1 fb~^ 
when the bunches are full. We will examine results achievable for m/j = 110 GeV 
and rrih = 130 GeV by employing only F = bb. 

The invariant amplitude squared for Higgs+7* + Z* exchange is given in the 
Appendix as a function of the relative angle ( between the transverse polarizations 
of the n~ and /i"*" and as a function of both 6 and 0, the angles describing the 
orientation of the b and b momenta in the final state. We first note that the 
interference term is never of importance for the cases we explore here. Secondly, 
we observe that \A4\ from Higgs exchange is independent of 6 and </>, whereas that 
for the 7* + Z* background does depend upon both 6 and (p. Our analysis proceeds 
as follows. 

We assume input values for (a, 6) = {a,b)/ {gm^/2mw) {a = a^ijb = 6^), de- 
noting them by ag and b^. We integrate over (p, but bin events in cos^^ (bin label 
j). As described, the spin precession means that we must also bin events in 6~ 
(remembering that 6^ will always have a known correlation with 6~ for a given 
spin-precession configuration); 9~ bins are labelled by i. We devote L/6, L/6 and 
2L/3 to the spin-precession configurations C = I, II, III described in the previ- 
ous section.0 The number of events (including both signal and background) for a 
given configuration choice C, given 6~ bin i and given cos^ bin j is denoted by 
N{C,i,j). The statistical error for each such bin is AN{C,i, j) = ^N{C, We 

"^In this paper, all R values will be quoted in absolute units and not in per cent. 
"^In all cases studied, if there are no a priori restrictions on the possible couplings, sensitivity 
to b/a is maximized by employing all the configurations. 
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then consider a different model characterized by a and b and compute 

^ j2 - N{C,i,j,do,bo))^ ^^^^ 

c,i,j AN'^{C,i,j,do,bo) 

The discussion of Section 2 shows that with appropriate binning and appropriate 
choices for our configurations we can effectively drop the C,i,j dependence of 
AN{C,i, j). The following discussion will implicitly rely on this fact. 

The accuracy with which d and b can be determined can be assessed by drawing 



contours of constant Ax^ in 6 = tan~^6/d, r = y a? + b'^ parameter space. The 
na error in 6 for example (which implicitly assumes that r is always adjusted 
for a given b/d so as to minimize the Ax^) is given by drawing the vertical line 
tangent to the Ax^ = contour. We will also include the Ax^ = 6.635 contour 



corresponding to 99% CL. The corresponding statistical errors for r = y a? + b"^ 
are obtained by the horizontal lines tangent to the contours. However, the true 
errors for r must incorporate the error coming from our imprecise knowledge of 
BR{h bb) and r*°*. 

Let us begin with the example of a SM Higgs boson (do = 1, 6o = 0) with 
mass TTih = 110 GeV. To illustrate the importance of polarization, we use as a 
reference point a total integrated luminosity at P = 0.2 (with full muon bunches) 
of L = 0.15 fb~^ (per year). We plot contours in (5, r) parameter space for some 
sample cases: (i) P = 0.2, 1 = 1 (i.e. nominal proton source intensity, >CpJ; (ii) 
P = 0.39, 1 = 1 (bunch merging at CpJ; (iii) P = 0.48, 1 = 2 (bunch merging 
at £ps = 2£pg); (iv) P = 0.45, / = 3 (full bunches at £ps = 3£pJ. The contours 

are presented in Fig. ^ It is useful to keep in mind that b/d = tan 5 is 1 for 
5 = 7r/4 ~ 0.785. For the cases (i)-(iv) above, the la, 2a, 99% CL, and 3(j 
limits on 5 are: (i) 0.94,. . .; (ii) 0.30,0.64,0.89,1.14; (iii) 0.20,0.41,0.53,0.64; and 
(iv) 0.15,0.32,0.42,0.50. The corresponding la and 99% CL limits on 6/d = tan 5 
are: (i) 1.36,. . .; (ii) 0.31,1.23; (iii) 0.20,0.58; (iv) 0.15,0.45. We see that even / = 1 
with bunch merging gives a reasonable la measurement; however, a good 99% CL 
limit on 6/d requires / > 2. 

Of course, as the Higgs mass increases, the event rates decline and the results 
worsen. This is illustrated in Fig. ^ which gives the same plots as Fig. ^ but for 
mh = 130 GeV. Above rrih = 130 GeV, the hb branching ratio begins to decline 
sharply and one should see a significant WW* final state rate, which in itself 
would be a signal of a significant CP-even component for the Higgs boson. A high 
discrimination check of the CP nature of the muonic coupling would be possible 
by including the WW* channel in the transverse polarization analysis. We will not 
pursue this here. 

Another perspective on our errors is to ask how well separated the b/d mea- 
surement for a SM Higgs boson is from the same measurement for an alternative 
Higgs boson of the general two-Higgs-doublet model (2HDM). To illustrate, con- 
sider nih = 110 GeV with do = &o = The Ax^ contours for this input model are 
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Figure 3: We give contours at Ax^ = 1,4,6.635,9 in the S,r parameter space 
assuming that the integrated luminosity at polarization P = 0.2 (with full muon 
bunches) is L = 0.15 fb^^. Four cases are compared: (i) P = 0.2, L = 0.15 fb~^ 
(which corresponds to / = 1); (ii) maintaining same proton source intensity, / = 
1, but merging the bunches, corresponding P"^(/ = 1) ~ 0.39, for which L — 
0.075 fb~^; (iii) increasing the proton source intensity by a factor of two, 1 = 2, 
while merging the bunches, corresponding to P™-{^I = 2) ~ 0.48, for which L = 
0.075 fb~^; (iv) / = 3, using just-full bunches, corresponding to P'^(/ = 3) ~ 0.45, 
for which L = 0.15 fb"\ We assume a SM Higgs boson (cio = 1, &o = 0) with 
ruh = 110 GeV and an overall efficiency (including 6-tagging efficiency) of 0.54. 
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Figure 4: Same as Fig. but for nih = 130 GeV. 
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Figure 5: Same as Fig. but for a mixed-CP higgs boson with Sq = 6o = 
and ruh = 110 GeV. 



18 



given in Fig. |^. The la contours are already non-overlapping with those for a pure 
CP-even Higgs (Fig. ^ for P = 0.39 and 7 = 1. However, the 99% CL contours 
have a small overlap even at / = 3 (P = 0.45). This again emphasizes the potential 
importance of larger-than-nominal source luminosity for such discrimination. 

Table 1: Event number pattern for different Higgs models as a function of (, 
assuming P^ = and P^ = P; see Eq. (|I]). 
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Let us now turn to a situation in which there are two degenerate Higgs bosons. 
This can arise in the 2HDM and also in the MSSM. The rates as a function of the 
relative angle between the transverse polarizations, (, will depend in detail upon 
the a and b values of the two Higgs bosons. As an example, suppose that we have a 
degenerate pair at ruh = 110 GeV, one of which is pure CP-even with a = 1/ ^/2, b = 
((5 = 0) and the other pure CP-odd with d = 0,b = 1/^2 {6 = n/2). Here, we 
have chosen the normalizations so that the total unpolarized (i.e. averaged over 
the four ( settings) production rate is the same as for a SM Higgs boson, with 
each of the two Higgs contributing equally to this rate. (We will also assume that 
the two Higgs bosons have the same bb branching ratio as a SM Higgs boson.) 
In this situation, Eq. (0) shows that the production rate will have no dependence 
on ( [cos(0 + C) + cos(7r + ^) = 0], whereas any single Higgs model will exhibit 
a distinct pattern as a function of (, as illustrated in Table |1|. Were we able to 
accumulate events at fixed (, the pattern of Ax^ for discriminating any two models 
from one another (assuming all have the same bb production rate) is apparent in 
the approximation where the ( dependence of the AN{() errors is neglected (as 
approximately appropriate given the dominance of the background contribution 
to AA^ and the weak dependence of the background on (); Ax^ would simply be 
proportional to the squares of the rate differences between two models summed over 
the four ( values. In the case where we compare the ^^independent rate predicted 
in a degenerate Higgs pair model to expectations for any given single Higgs boson, 
it is apparent from the table that the Ax^ between the degenerate model and any 
single Higgs model is independent of the latter. (This is true for arbitrary 6 for a 
single Higgs since cos2(25) + cos2(25 + tt/2) + cos2(25 + tt) + cos2(2(5 + 37r/2) = 2.) 

However, because of spin precession, we must actually employ Eq. (^. Since 
sin 25 = for both Higgs bosons while cos 25 = 1 for the CP-even and cos 25 = — 1 
for the CP-odd Higgs boson, after summing over both, the cross section will be 
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the same as that obtained if we set cos 25 = as well as sin 25 = 0. To compare 
to the simplified discussion of the previous paragraph, let us recall that for our 
choices of configurations and binning we can approximately neglect the C,i,j (i.e 
configuration, 6~ and 6) dependence of the AN{C,i, j) denominators in Eq. (|2l|) . 
Then, for the L/6, L/6, 2L/3 luminosity weightings for configurations C = I, C = 
II, C = III, respectively, the effective sensitivity for discriminating between the 
[cos 25 = sin 25 = 0] — equivalent situation for the degenerate pair and the results 
expected for any single Higgs boson characterized by angle 5 will be proportional 
to |<S2 + ~ ^ _ ip4 |^pQg2 2S + sin^ 25) = [see Eqs. (|TB|) and 

(^Ol) ], i.e. again independent of the 5 value for any single Higgs boson to which one 
might compare. One finds that this approximation is actually quite good. Taking 
L = 0.15 fb"^ for J = 1 and P = 0.2 and assuming a SM-like bb production rate, we 
find Ax^ = (i) 0.38; (ii) 2.8; (iii) 6.4; (iv) 9.8, for the four (/, P) scenarios defined 
earlier, essentially independent of the type of single Higgs boson exchange to which 
one compares. Note that we get > 99% CL exclusion of a single Higgs model only 
for J > 2 if nature chooses a degenerate pair. In any case, the ( dependence is 
key to separating a degenerate CP-even plus CP-odd pair of Higgs bosons from a 
single Higgs boson of any given CP nature. 

In the MSSM, in the absence of CP violation the most likely situation in which 
we will encounter a highly degenerate pair of Higgs bosons that cannot be easily 
separated by scanning in ^/s is in the limit of large m^o and large tan/5. The 
increasing degeneracy with increasing tan f3 is illustrated for rriAQ = 400 GeV in 
Fig. assuming squark masses of 1 TeV and no squark mixing and a beam energy 
spread of R = 0.001. Since the total widths of the if" and are substantial 
(> 1 GeV) for the m^o and tan [3 values being considered, it is not guaranteed 
that we will be able to separate the peaks. The figure shows that we are able to 
observe two separate peaks (the peak being at lower mass than the peak) 
for moderate tan/? < 6. But, for higher tan/5 values the peaks begin to merge; 
for tan/3 > 8, \mHO — m^o| < 1 GeV and one sees only a single merged peak. The 
picture changes if squark mixing is substantial; for instance, for m^o = 300 GeV, 
squark masses of 1 TeV and large squark mixing {At = Af, = 3 TeV), the and 
A° peaks actually cross at tan /? ~ 5. 

To explore the various possible scenarios, we begin by considering our ability 
to discriminate between an exactly degenerate pair of CP-even and CP-odd Higgs 
bosons vs. a single Higgs boson (of a given type) as a function of tan p. We note 
that the present situation is significantly more favorable than that discussed above 
with a light degenerate pair with SM-like widths, branching ratios and production 
rate. Since the widths of the and A^ are significant, one can operate the 
muon collider with the natural beam energy resolution of order R = 0.001 and 
still have < T, which maximizes the Higgs production rate.0 At such R and 

''^Of course, one could possibly separate the two Higgs bosons by employing R < 0.0001, but 
it is unlikely that the machine would be operated in this way at this higher energy unless a Higgs 
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Figure 6: We plot bb (solid) and ti (dashed) event rates for total integrated lumi- 
nosity oi L — 7 pb~^ coming from h'^/jl" — > + A° as a function of ^/s, assuming 
m^o = 400 GeV. Each window is for the specific tan P value noted. These event 
rates are to be multiplied by a factor of 1000 for the expected yearly integrated 
luminosity of 7 fb~^. We employ squark masses of 1 TcV and no squark mixing. 
Supersymmetric decay channels are assumed to be closed. 
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for \/s ~ 300 GeV (400 GeV), the nominal yearly integrated luminosity for full 
bunches at P = 0.2 is estimated at L ~ 2 fb~^ (7 fb~^), i.e. more than a factor of 
ten larger than the R = 3 x 10~^ value. Second, the bb branching ratio for large 
m^o and large tan (3 is inevitably of order 88%, the only significant competitor 
being r+r^. Even for moderate tan/5 and for Higgs masses above ti threshold, the 
bb branching ratio remains substantial. Even more importantly, at large ttiao the 

IjT couplings of the and are enhanced relative to SM strength by a factor 
of tan/9 so that the fi~^fj^~ branching ratio asymptotes to a constant. 

As shown in Figure ^, the and need not give exactly the same bb event 
rate (due to differences in BR{^^ ^^)BR{hb)). This is true even if parameters are 
chosen so that they are exactly degenerate in mass at some given tan /3 value. On 
the other hand, it is also typically possible to choose parameters so that they do 
have exactly the same rate. We consider this last possibility first. (As already 
discussed, when the and A^ have exactly the same hb rate, there is no depen- 
dence of the summed hb rate onC,. As a result, even after including spin precession, 
Ax^ will be independent of the CP nature of any single Higgs boson to which we 
compare if we sum over the three configurations C = I, II, III with luminosity 
weighting L/6, L/6, 2L/3.) To simulate this situation, we choose a value of tan/? 
and a value of the A^ mass. We then compute the and A^ event rates (assuming 
stop masses of 1 TeV and no squark mixing) and reset the and A^ event rates 
so that both are equal to the average of the originally computed and A^ rates. 
We consider m^o values of 300 GeV and 400 GeV. In the latter case, the ti decay 
channel is open and has significant branching ratio at lower tan (3 values. This, 
along with the decreased coupling at lower tan (3, is included in computing 

the bb rate. We make no use of the ti channel; its inclusion would, of course, in- 
crease our discrimination power, especially if the final state correlations that can 
be probed there are employed. The resulting discrimination power is generally 
more than adequate without its inclusion. 

Figs. 1^ and ^ give the results for m^o = 300 GeV (400 GeV), respectively. We 
plot Ax^ obtained for the four polarization/proton-source-intensity options (i)- 
(iv) delineated earlier, but using the higher nominal luminosities stated above: 
(i) P = 0.2, and the / = 1 values of L = 3(10.5) fb"^ (ii) P = 0.39, and the 
1 = 1 merged-bunch L = 1.5 (5.25) fb~^ values; (iii) P = 0.48 and / = 2, yielding 
merged-bunch L = 1.5 (5.25) fb^^; (iv) P = 0.45 and / = 3, without bunch 
merging, yielding L = 3 (10.5) fb"^. We emphasize that options (i) and (ii) require 
no over-design of the proton source. The Ax^ plots show that good discrimination 
is obtained even for option (i) once tan/? > 10. Option (ii) would be needed for 
good discrimination if tan/3 ~ 5. 

We next consider the discrimination power achieved if we continue to enforce 

resonance peak is seen and there is already some evidence, through the techniques considered 
here, that there are actually two overlapping resonances with different CP properties. Even then, 
there will remain the possibility of very close or even exact degeneracy. 
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Figure 7: In the upper left window, we plot the bb event rates in the MSSM for 
the H^^ and (the rate is the larger of the two) as a function of tan/3 for 
m,A0 — 300 GeV, assuming squark masses of 1 TeV, no squark mixing and inte- 
grated luminosity of L = 2 fb"''^. Also shown is the (relatively small) background 
rate. In the remaining windows we plot A^^, after including precession and in- 
creasing L to L = 3 fb^^, as a function of tan/? for three different cases in which 
we forcibly lower uiho to 300 GeV (for exact degeneracy). (1) We adjust the 
and A° event rates so that each is exactly equal to the average of the ^4° and 
rates as predicted by the MSSM, and compute Ax^ for H° + A^ vs. a single 
Higgs resonance (of any type) with the same total event rate, employing only the 
bb channel. (2) We use the actual event rate and compute Ax^ for + A^ vs. 
a single CP-even resonance yielding the same bb event rate. (3) As in (2), but vs. 
a single CP-odd resonance. In cases (l)-(3), we give results as a function of tan/3 
for the four polarization-luminosity situations (i)-(iv) (as labelled on the curves) 
described in the text. no 




Figure 8: The same as Fig. ^ except for uiao = mno = 400 GeV, and using 
L = 7 fb~^ for the rate window increased to L = 10.5 fb~^ for the Ax^ windows (in 
order to account for the inefficiency associated with spin precession). The upper- 
left window gives the and and background tt event rates as dashed curves. 
In the Ax^ analysis, only the bb final state (after accounting for the depletion from 
H^,A° ^ tt) is employed. 
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Figure 9: We plot contours of Ax^ = 1,4,6.635,9 in the A, ^2 parameter space, 
assuming m^o = 400 GeV, tan/3 = 7, squark masses of ITeV and no squark 
mixing. The different windows give results for the different proton source intensity 
and bunch merging options described in the text. 



25 



<] 




5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10 

tan^ 



Figure 10: We plot A^;^ as a function of tan (3 for discriminating between the 

mixture predicted by the MSSM model (as specified in the text) at = rnuo 
compared to a single CP-even or CP-odd Higgs boson yielding exactly the same 
bb event rate. Results are presented for the polarization/proton-source-intensity 
options (i)-(iv) described in the text. 
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exact degeneracy by lowering uiho to m^o, but employ the actually predicted A" 
and event rates. For m^o = 300 GeV, the event rate for bb is lower than the 
A° rate, and the worst (best) discrimination is achieved relative to a purely CP- 
odd (purely CP-even) Higgs boson (again assumed to have exactly the same total 
bb event rate). The resulting Ax^ values for options (i)-(iv) are given in the two 
bottom windows of Fig. |^. Even in the worst case, good discrimination power is 
achieved for tan /9 > 5 by using option (ii). Exactly the reverse situation arises for 
rriAo = 400 GeV. As seen in Fig. ^ the rate in the bb channel is smaller than for 
the H^. Discrimination against a purely CP-even Higgs will be the most difficult. 
In fact, we see that option (ii) will not provide clear discrimination against a purely 
CP-even Higgs if tan P <7. 

Finally, let us consider the case where the and are not exactly degenerate, 
but are only nearly so. To be specific we adopt the MSSM predictions of Fig. ||. 
As noted earlier, if the collider is operated with R = 0.001 beam energy spread 
so as to maximize luminosity, for tan (3 > 8 one observes a single broad peak in a 
scan and we will not know that there are two Higgs bosons present. The apparent 
peak of the resonance shape will tend to coincide with the mass of the Higgs boson 
which yields the higher bb rate (which Higgs this is depends on the mass, as we 
have seen) . If we center on this apparent peak, the contribution to the bb rate from 
the weaker Higgs will be further reduced because ^/s is somewhat on the wing of 
its resonance shape. This further decreases the Ax^ discrimination power relative 
to a purely CP-odd (CP-even) Higgs boson if the dominant Higgs is CP-odd (CP- 
even). Even if we see two separate peaks, we will wish to experimentally determine 
which is the and which is the A^. _ 

To illustrate, let us define = \/ N^o + Njro and A = tan""*^ where N^o 

{Nho) is the number of A^ {H^) events at the chosen ^/s for a possible model divided 
by {N^o + N^o)]^^ggj^i- In Fig. ||, we present Ax^ contours in A,Q parameter 
space for the input MSSM model (see Fig. 1) specified by m^o = 400, tan (3 = 7, 
squark masses of 1 TeV and no squark mixing. To determine that there are 
events under the peak at 99% CL requires the 1 = 2 proton source intensity 
option (iii). A global overview of our ability to discriminate the MSSM input 
model from a single purely CP-even (the worst case) or single purely CP-odd (the 
best case) Higgs boson yielding exactly the same bb event rate is provided by 
Fig. |lO[ There, we plot Ax"^ for these two discriminations as a function of tan /3 for 
polarization/proton-source- intensity options (i)-(iv). We find that discrimination 
against a purely CP-odd Higgs boson is excellent even for tan/3 ~ 5, whereas 
tan/3 > 9 [tan/5 > 7.3] is required for 99% CL discrimination against a single 
CP-even Higgs boson using option (i) [(ii)]. Of course, if squark mixing is large, 
the degree of degeneracy between the and A^ can be such that the peaks 
will merge even when tan/? ~ 5. In this case, good discrimination power would 
typically require enhanced proton source intensity. 
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5 Summary and Conclusions 



The most natural polarization for the muon bunches at a muon collider is P ~ 0.2. 
For this polarization, the collider luminosity will be maximal. In initial exploration 

of the Higgs or other narrow resonance, one will choose the polarizations to lie 
in the horizontal plane of the storage ring. As the polarizations rotate it will 
then be possible to perform a precise measurement of the beam energies and their 
Gaussian widths and of the degree of polarization itself. Further, by choosing the 
/x"*" polarization to be 90 degrees out of phase with the /i" polarization, the effect 
of the rotating polarization will cancel out when averaging over the 1000 turns 
during which the typical bunches are stored; i.e. the turn-averaged cross section 
will be identical to the polarization-averaged cross section. However, once the basic 
properties of the Higgs boson or resonance are known (total width, branching ratios, 
etc.) the next important goal will be to determine its CP properties as codified in 
the relative strength of its scalar and pseudoscalar couplings to fermions. The muon 
collider provides a perfect opportunity for determining these relative strengths for 
the muon itself, but only if one has retained the ability to reconfigure the collider 
so as to run with high polarizations for the muon bunches and to have flexibility in 
the orientation of these polarizations at the time the bunches are inserted into the 
storage ring. In most muon collider designs, high polarization can only be achieved 
by making strong momentum cuts on the muons, which, unless the proton source 
has 'spare' luminosity relative to bunch saturation limits for P ~ 0.2, will result in 
some loss of coUider luminosity. However, if the goal is to measure the CP-odd/CP- 
even coupling ratio, this loss of luminosity is more than compensated by increased 
sensitivity. Our goal in this paper has been to develop efficient techniques and 
polarization configurations for determining the CP-odd/CP-even coupling ratio 
and to quantify the accuracy with which this ratio could be extracted from the 
data. 

We have found that one very effective technique is to accumulate events for 
three carefully chosen polarization configurations (as defined by the polarization 
of the yU"*" and fi~ bunches at the time they enter the storage ring). To achieve 
these three configurations, one will need appropriate solenoids and/or small rings 
for manipulating the polarizations prior to injection into the storage ring. Once 
the bunches are in the storage ring, further manipulation would destroy our ability 
to measure the bunch energies to the 1 part in 10^ level needed (at least for a 
narrow Higgs boson). We have demonstrated that by selecting high-polarization 
muons and performing bunch merging, meaningful constraints on the CP nature 
of the muonic couplings of a Higgs boson are possible even if it is light and narrow 
and even if the proton source has only the nominal luminosity required to saturate 
the bunch limits in the storage ring at relatively low polarization. However, we 
have seen that extra proton source luminosity (sufficient to saturate the bunch 
limits of the storage ring when muon selections leading to large polarization are 
employed) may be needed to achieve a high level of certainty regarding the CP 
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nature of the muonic couplings of such a Higgs boson. We have also shown that 
distinguishing a highly degenerate pair of CP-even and CP-odd Higgs bosons from 
a single Higgs with definite CP nature will generally be possible, especially if this 
is the degenerate H^-A^ pair of the MSSM at high tan j3 and large 171^0 for which 
event rates are high and high-luminosity-running at i? = 0.001 would suffice. 

Overall, while the machine capabilities required to make a good determina- 
tion of the CP nature of the muonic couplings of a Higgs boson do not come 
without a price, it could well happen that the true nature of an observed Higgs 
resonance would be obscure without the measurements considered here. Further, 
Higgs/resonance CP studies provide but one example of how unique sensitivity 
to new physics can result if we can take advantage of the slow precession of the 
muon bunch polarizations by manipulating their orientation and relative phases 
and recording events on a turn-by-turn basis. Thus, we encourage the muon col- 
lider designers to retain the flexibilities needed to insert devices for rotating the 
spins of the beams to a variety of initial (i.e. prior to insertion into the storage 
ring) configurations of longitudinal and transverse polarization and to consider 
seriously the possibility of 'over-designing' the proton source relative to storage 
ring bunch saturation limits associated with low polarization. Or, perhaps it will 
prove possible to employ a very high gradient in the initial phase-rotation stage of 
pion/muon capture JlT], [l^, thereby maximizing the luminosity available after 
the selection cuts required for high polarization. 
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6 Appendix 

In this appendix, we present the analytic formulae for the signal and background 
cross sections in an arbitrary fermionic final state. Fermionic couplings take the 
form i 7^6/ (photon), i'-i^{cf + dp^) (Z), and i{a f + i'y 5b f) (Higgs). Here, e/ = eQf, 
= 2^^iTi - 2Qfsm'^9w), dj = -j^^Ti, and, for a SM Higgs boson, 
a/ = and bj = 0. In terms of these quantities, the decay width of the Higgs 
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boson to // is given by 

nh - //) = ^f3f m, {a}(3j + b)) , (22) 



57r 



where jSj = \jl — 

We employ the center of mass system with -y/s = E'cm and define C = C^~C~ as 
the angle of the /x"*" transverse polarization relative to that of the The explicit 
expressions for the 4-momenta and spin vectors of the /i"*" and /i~ in the laboratory 
frame are: 

p^- = ^(l,0,0,/3) 
p^+ = ^(1,0,0,-/?) 

s^- = PE^i(3, 0, 0, 1) + (0, cos C~, sin C", 0) 
s,+ = P+7(/5, 0, 0, -1) + P+(0, cos C+, sin C+, 0) 
/s 

Pf = -^{1, Pf sinO cos (j), /3f sm6 sin (f), /3f cos6) 

Pf = "^(l; ~Pf sin 6' COS 0, —f3f sin 6' sin 0, — /5/ cos6') , 

(23) 

where (j) and 6 are the standard angles defined using polar coordinates for the final 
fermion. The forms for s^- and above are those which make the separation 
between longitudinal and transverse polarization most evident. The conversion 
between these forms and those given earlier in Eqs. (|^) and appropriate to 
making the precession physics most transparent can be accomplished by the fol- 
lowing mappings: 

P| cos = P^, -P| sin = COS , Py = Pt sin • (24) 

We now give expressions for the invariant matrix element squared, |A^|^, for 
H^jj,^ ff, summed over final spins and averaged over initial spins. These will 
be given in terms of P^, P-^ and C^. Eq. (|2^) can be used to convert to the P^, 
Py and 6^ precession variables. We divide \Aif into three pieces: the absolute 
square of the Higgs diagram, the interference between the Higgs diagram and the 
7* + Z* exchange background diagrams, and the absolute square of the 7* + Z* 
background diagrams. The symbols 3? and 55 denote the real and imaginary parts. 
We give the results for m^/ ^/s —>■ 0, an excellent approximation for any reasonable 
Higgs mass. Defining Ux = — m'j^ + iTxmx] we have the following. 
Absolute square of Higgs exchange: 

(aj/3j + bj) \U,\' s' [(aj + fej) (1 + Pf P+) + Pf P+ { (aj - 6j) cos C - 2 a, b, sin (} 
= {a'f/3j + b)) |H,|^ (al + 6j) [(1 + P^P^) + P^ P} cos(C + 25)] (25) 
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Higgs - 7* + Z* Interference: 



1 /2 

- sin(C+ + 5-<P)P^ eje, 53(n^n*) + cj (c^ + d, P^) %^z^l] 



— sin(C — 5 



T 



ej e, %U,ni) + Cf {c, - d, Pt) 53(n^n* 



- cos(c+ + 6-<p)p^ [cfe^ p^ 5R(n^n*) + Cf (c^ p^ + d^) K(nzn*) 
+ cos(c- -5-<p)p^ [eje^ p+ 3fJ(n^n*) + cj {c, p+ - d,) 3?(n^n*) 

Absolute square of 7* + Z*: 



(26) 



-8 4 



2 ITT |2 

rUf \LLz\ s 



{cl + dl) (1 - P^Pl) +2c,d, (P, - Pi 
+2 sin(C - 20) sin^ 9 P^ P^Cf d^ e/ "^{Tiz) n^, /5j 
+ cos(C-20) sin2^^P^P+|(4 + 4) |nz|' + e5ejn2 + 2e;e^c;c^n^3ft(n, 

+ (2 - /3f sin^ ^) I (4 + 4) [(cj + 4) (1 - P^Pt) +2c,d, {P, - Pi 



+2 Cf Cf 



1 - PlPE 



d, {PE - P^)] mz) + e} el (1 - p^-p+) n^l 



+4 cos^^ rf/jc/ [2c^f/^ (1 - P^-P+) + (cj + 4) (P^- - P^ 
+e; [d^ (1 - P^-P+) + c, {PE - Pi:)] n, mz)]Pf s' 

In terms of |A1p, the cross section as a function of s, 6' and is given by 
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d cos 6' ( 



6471^5 



Pf\M\' 



(28) 



For example, if we spin average the Higgs portion of \Ai\'^ by combining C = 
and ( = Tc (so that the average of cos^ is zero), integrate over dcosO and d(j), 
and use Eq. (^), then we obtain the standard spin-averaged total cross section: 
(j(s) = 47rT{h — > n~^n~)T{h — > //) \Ilhf- When this latter form is convoluted with 
a Gaussian distribution in y/s, one obtains the result of Eq. (^. 
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